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ABSTRACT 

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if 
accretion onto Supermassive Black Holes (SMBHs) and star formation are explicitly linked up to z ~ 1. 
We identify 152 galaxies that harbor AGN, based on their X-ray emission (Lo.5— 10 keV 

> 10 42 erg s" 1 ) 

detected by X M M — Newton observations of 7543 galaxies (i aC s < 22.5). Star formation rates (SFRs), 
including those weighted by stellar mass, of a subsample are determined using the [OII]A3727 emission- 
line luminosity, corrected for an AGN contribution based on the observed [OIIIJA5007 strength or that 
inferred by their hard (2-10 keV) X-ray luminosity. We find that an overwhelming majority of AGN 
host galaxies have significant levels of star formation with a distribution spanning ~ 1 — 100 M Q 
yr _1 ; their average SFR is higher than that of galaxies with equivalent stellar mass (M» > 4 x 10 10 
Mq). The close association between AGN activity and star formation is further substantiated by an 
increase in the fraction of galaxies hosting AGN with the youthfulness of their stars as indicated by 
the rest-frame color (U-V) and spectral index D n (4000); we demonstrate that a mass-selected sample 
is required to alleviate an artifical peak in the AGN fraction falling in the transition region due to the 
fact that many 'blue cloud' galaxies have low mass-to- light ratios in luminosity- limited samples. We 
also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the 
overall galaxy population and naturally explains the low SFRs in AGNs (z < 0.3) from the SDSS. We 
conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large 
reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR 
is shown to be weak although the average ratio (~ 10 -2 ) is constant with redshift, effectively shifting 
the evidence for a co-evolution scenario in a statistical manner to smaller physical scales (i.e., within 
the same galaxies). The order-of-magnitude increase in this ratio compared to the locally measured 
value of Mbh /Mbuige, is consistent with an AGN lifetime substantially shorter than that of star 
formation. Our findings illustrate an intermittent scenario with underlying complexities regarding 
fueling over vastly different physical (and temporal) scales yet to be firmly determined. 
Subject headings: quasars: general, galaxies: active, galaxies: evolution, X-rays: galaxies, quasars: 
emission lines 
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IShankar et ail l2007t ISilverman et~ail I2008H) suggests 
a common mechanism that regulates their growth . 
Major mergers of gala xies (iMihos fe Hernquistl Il996f) 
and secular evolution (|Kormendv fc Kennic utt 2004; 
iHopkins fc Hernquistl 120061 ) can both potentially funnel 
gas to the nuclear region that can then power concurrent 
star formation and accretion onto supermassive black 
holes. Fully consistent with this scheme, young stellar 
populations are known to be prev alent within the bulges 
of nearby Seyfert galaxies (e . g.. iTerlevich et al.1 IT990I ; 
iGonzalez Delgado et al.l 120011: IWhittle fe NelsonT 12003ft 
and l uminous quasars ( Jahnke et al.ll2004al : iLetawe et alJ 
120071 ). Enhanced levels of ongoing star formation 
are also evident in nearby Seyfert 2s (|Gu et alj l2006f ) 
and the hosts of some mid-infrared selected quasars 
(|Lacv et al.l 120071 ) although are not high enough to 
put them in the same class as the Ultraluminous In- 
frared Galaxies (|Sturm et al.l [2006) - On the other 
hand, the lack of star formation based on [Oil] emis- 
sion in_the hosts of a significant sample of PG q uasars 
(flol I2005T ) and type 1 AGN ijKim et alj 120061 ) from 
the Sloan Digital Sky Survey (SDSS) may indicate an 
elevated role for AGN in suppressing s t ar formation 
(iGranato et al.ll2004t ISpringel et aLll2005t ICroton et all 
120061: ICavaliere fc Mencil 120071: IHopkins et all 12008a) . 
However, there do appear to be high levels of star 
formation in the hosts of PG qua sars when con- 
sider i ng their far- i nfrare d emission ([Schweitzer et all 
120061: iNetzer et all l2007[ ). and type 2 quasars from 
the SPSS (0.3 < z < 0.8; iZakamska et~ail l200l 
iKim et ail 120061 ) that may hint at an underlying evo- 
lutionary sequence (i.e., a modification to the uni- 
fied model with a correlation between star formation 
and nuclear obscuration) or the systematic increase 
in the global star formation history of galaxies (e.g., 
Lilly et al. 1995; Hopkins & Bcacom 2006; Nocskc ct al. 



20071 iTresse et al.H ~2007: Zhe ng et al.ll2007t ). These seem- 



ingly discrepant results highlight the complexity in de- 
termining if a casual relationship between star formation 
and AGN activity exists. This is most likely due to the 
challenges in measuring the host galaxy properties un- 
der the glare of a bright AGN, widely varying selection 
methods, and inadequate control samples. 

The ability to study the host galaxies of AGN se- 
lected from a large parent sample of galaxies has dra- 
matically improved in the last few years. The SDSS 
has generated an unprecedented sample of galaxies to 
select larg e numbers of low redshift (z < 0.3) narrow- 
line AGN dKauffmann et al.ll2003bl: iKewley et "ail 120061: 
Stasihska et all 120061 : iHeckman et all 120061 : I Wild et al I 



2007h to cleanly study their host properties. Pow- 
erful AGN have been clearly shown to reside in 
bulge-dominated galaxies with young stellar populations 
(jKauffmann et al.ll2003bf ). In conjunction with GALEX 
observations, AGN have optical- ultr aviolet colors plac- 
ing them in an in t ermediate region (iMartin et al.l 120071 : 
ISalim et al.l 120071 : ISchawinski et al.l I2007T ) between the 
red and blue galaxy populations thus s tirring interest 
in th eir relation to tra nsitional galaxies (IHopkins et al.l 
I2007fh Remarkably, iKauffmann et all (|2007j) find a 
strong correlation between the mass accretion rate onto 
SMBHs and the presence of young stars in the bulge 
that is typically accompanied by a star-forming outer 
disk that may be supplying the nuclear fuel. It is imper- 



ative to extend such studies to epochs closer to the peak 
[z ~ 1) of the cosmic AGN and star formation history 
since AGNs may have an impact on the gas content of 
galaxies and their subsequent demise. 

X-ray emission is a unique probe of AGN activity 
due to its ability to penetrate a moderately obscuring 
medium (Njj < 10 24 cm -2 ), the lack of confusion with 
stellar sources, and because it allows a more direct probe 
of the mass accretion rate (see lBrandt fc Hasingcr 12005L 
for a review). Selection of narrow line AGN (i.e., ob- 
scured; type 2) b y 
Szokolv et all 12004 



X-rays (e.g., 


Bargcr ct al. 2003; 


Silverman et al. 


20051: iBrusa et all 



20071 ) can compensate for the deficiencies of optical se- 
lection based on emission line ratios at z > 0.3. X- 
ray observations with both Chandra and XMM- Aeu>fon 
are being exploited to cleanly study the host galaxies of 
X-ray selected AGN dGrogin et all 120051: iNandra et all 
20071 iPierce et all 120071: iGabor et alj l2008fh Recently 



Silverman et al. (2008a) show that the bulge-dominated 



hosts of X-ray selected AGN in the Extended Chandra 
Deep Field - South have rest-frame colors that are bluer 
with increasing redshift, possibly related to the star for- 
mation history of galaxies. In a related study, the mean 
star formation rate based on broad-band photometry of 
58 X-ray selected AGN (z ~ 0.5-1.4) in the 1 Msec CDF- 
S has been shown to be similar to IRAC-selected galax- 
ies in an equivalent mass regime (jAlonso-Herrero et al.1 
120071 ). To date, few studies using a significant sample 
of AGN selected from a large parent sample of galaxies 
with optical spectroscopy have been undertaken at these 
higher redshifts (z > 0.3). 

To do so, we utilize the ri ch multiwavelength observa- 
tions of the COSMOS field (|Scoville et al.ll2007l ) to carry 
out such a study at 0.5 < z < 1.0. The COSMOS sur- 
vey is roughly a 2 square degree region of the sky se- 
lected to be accessible from all major observatories both 
from the ground (e.g., Subaru, VLT) and space (e.g., 
HST, Spitze r, XMM-Newton , Chandra). The zCOS- 
MOS survey ([Lilly et al.ll2007l ) targets objects for opti- 
cal spectroscopy with the VLT in two separate observing 
programs. A bright sample (i < 22.5) is observed with a 
red grism to provide a wavelength coverage of 5500 — 9500 
A ideal for identifying galaxies (£*) up to z ~ 1.2. A 
deep program, not utilized in the present study, targets 
faint galaxies (B < 25), selected to be in the redshift 
range 1.5 < z < 2.5, using a blue grism with a wavelength 
coverage of 3600 < A < 6700 A. Here, we select a sam- 
ple of galaxies based on their stellar mass with reliable 
spectroscopic redshifts from the zCOSMOS bright pro- 
gram. Those that host AGN are identifie d by their X-ray 
emission as detected by XM M-Newton (Hasingc r et al.l 
120071: ICappelluti et aTll2007f ). We measure the strength 
of emission lines (i.e., [OII]A3727, [QIII1 A5007), using our 
autom ated pipeline (" platefit_vimos" ; lLamareille et al.1 
2008) to determine star formation rates for the entire 
galaxy sample including those hosting AGN. An addi- 
tional spectral indicator (Z?„4000) enables us to discern 
the age of the stellar populations on longer timescales. 
X-ray luminosities of the AGN give us a handle on their 
bolometric output, less affected by obscuration, to infer 
mass accretion rates and determine any trends with star 
formation rate and redshift. Finally, we p oint out that a 
companion paper (Silverma n et al.ll2008cD based on the 



Star formation in AGN host galaxies 



3 



zCOSMOS 10k catalog expands on the current study by 
investigating the environmental impact on AGN activity 
and their host galaxy properties. 



Throughout this work, we assume Hq 
Mpc" 1 , n A = 0.75, and fi M = 0.25. 



70 km s 



2. DATA AND DERIVED BROAD-BAND 
PROPERTIES 

2.1. Parent galaxy sample and AGN identification 

We use the zCOSMOS 10k spectroscopic 'bright' cat- 
alog to construct a well-defined sample of galaxies in- 
cluding those hosting X-ray selected AGN up to z ~ 1. 
Specifically, we identify 7543 galaxies with a selection 
magnitude %acs < 22.5 and high quality spectra hence 
reliable redshifts up to z = 1.02. A galaxy is included in 
our sample if the redshift has a quality flag 2.0 or higher 
that amounts to a confidence of ~ 99% for the overall 
sample. The redshift success rate is not strongly depen- 
dent on galaxy color over 0.5 < z < 1.0 given the qual- 
ity of the spectra and presence of strong features (i.e., 
4000 A break; Ca H+K, [Oil]). Full details on data ac- 
quisition, reduction, redshi ft measurements and q uality 
assurance can be found in iLillv et all (120071.120081). 

We use the catalog of X-ray sources (|Cappelluti et al.l 
120071 ) generated from the uniform XMM.- Newton 
coverage (~ 50 ks depth) of the COSMOS field 
(|Hasinger et al.ll2007[) to identify galaxies within our par- 
ent sample that harbor AGN. The X-ray catalog includes 
1848 point-like sources detected above a given thresh- 
old using a maximum likelihood method in either the 
soft (0.5-2 keV), hard (2-10 keV) or ultra-hard (5-10 
keV) bands down to a limiting flux of 5xl0~ 16 , 2xl0~ 15 
and 5xl0 -15 erg cm~ 2 s _1 in the respective band. The 
adopted threshold (Likelihood > 10) corresponds to a 
probability ~ 4.5 x 10~ 5 that a source is a background 
fluctuation. An additional 26 faint XMM sources coinci- 
dent with diffuse emission (Finoguenov et al. in prepa- 
ration) were excluded. Both the soft (0.5-2.0 keV) and 
hard band (2.0-10.0 keV) detections comprise our AGN 
sample in order to include the low-to-moderate luminos- 
ity and/or absorbed sources. The higher sensitivity of 
the XMM- Newton observations in the soft band enables 
us to probe lower luminosity AGN (log Lx ~ 43) with 
0-8 z < 1.0 not sufficiently represented in the hard- 
band catalog. 

Optical and near-infrared counterparts to the XMM- 
Newton X-ray sour ces are found usi ng a maximum- 
likelihood method (jBrusa et alj I2007D . Most X-ray 
sources (84%) have reliable counterparts while the re- 
maining are uncertain due to the presence of multiple 
objects within the X-ray error box that have similar 
probababilities of being the true counterpart. We use 
the Chandra observations (Elvis et al. in preparation; 
Civano et al. in preparation) that cover the central 1 
square degree to further refine our identifications. We 
elect to include counterparts with lower confidence given 
that ~ 50% of them are likely to be correct based on 
the overlap between XMM-Newton and Chandra iden- 
tifications (M. Brusa, private communication). All re- 
sults presented here are confirmed based on the slightly 
smaller catalog of highly reliable counterparts. 

The zCOSMOS 'Bright' program provides optical spec- 
tra for 357 of the 1093 optical counterparts to X-ray 
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Fig. 1. — X-ray luminosity (0.5-10 keV) versus redshift for 152 
AGN (large symbols). Filled circles are AGN primarily charac- 
terized by narrow emission lines while open symbols denote those 
with at least one broad optical emission line falling within the ob- 
served spectral range. The upper limit to the broad-band X-ray 
luminosity for each zCOSMOS galaxy is shown by the small grey 
dots. The solid horizontal line is the minimum X-ray luminosity 
that we enforce to measure the fraction of galaxies hosting AGN 
to ensure a significant underlying galaxy population. The dotted 
line shows our chosen upper limit to avoid AGN contamination. 



sources having i < 22.5. Reliable spectroscopic redshifts 
are available for 90% of which 164 have redshifts between 
< z < 1.02. Higher redshift AGN up to z ~ 4 have 
been identified by the zCOSMOS 'Bright' program but 
are not the focus of this study. Since a high fraction of the 
zCOSMOS galaxies with associated X-ray emission have 
detections in the hard 2-10 keV band (84%), we do not 
expect the inclusion of soft-selected sources to induce any 
significant bias. A subset (54%) of the X-ray sources has 
been designated as 'Compulsory' targets when designing 
the VIMOS masks which essentially lessens their ran- 
domness and results in a 71.9% sampling rate compared 
to the current rate of 29.8% for galaxies. We set priority 
to these sources since they were not include d in the spec- 
trosc opic followup program using Magellan (| Trump et al.l 
I2007I) . Since the number of AGNs is not overwhelmingly 
large, we chose to not select a smaller pseudo-random 
sample but rather incorporate correction factors applied 
to derived measurements where applicable. 

Moderate-luminosity AGN (L x ~ 10 43 erg s -1 ), de- 
tected in medium-to-deep X-ray observations, have been 
shown to be ideal laboratories for the study of their 
host galaxies due to their low optical brightness which 
in many cases is due to their nuclear obscuration (e.g., 
iTozzi et al.ll2006HMainieri et al.ll2007l ). Our AGNs pre- 
dominately have luminosities below that of typical QSOs 
(Lx > 10 44 erg s _1 ) but are the dominant contri butor 
to the Cosmic X-ray Background (jGilli et al.ll2007l ). In 
specific cases, we further select AGN based on their X- 
ray luminosity (42.0 < log Lq. 5-10.0 keV < 43.7), as done 
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in lSilverman et aT] (|2008a|) . to isolate a sample for which 
we can cleanly determine their host propert ies (i.e., stel- 
lar masses, rest-frame U - V, L>„4000). iGabor et all 
(2008) have demonstrated that the optical emission from 
a nearly equivalent AGN sample in COSMOS is primarily 
attributed to their host galaxies; AGN dominated galax- 
ies are preferentially at z > 0.8 due to the fact that 
their X-ray luminosities are typically above our cutoff 
(Lx ~ 10 43 - 7 erg s _1 ). We conclude that the derived 
properties of the host galaxies of AGN in zCOSMOS are 
likely to be reliable. Where feasible, we include higher 
luminosity AGN (QSOs; 16 with L .5-io kcV > 10 44 
erg s _1 ), since we can measure their stellar properties 
(i.e., [Oil] strength) even under the glare of an intrinsi- 
cally, bright AGN (see Section [3] for details). For these 
cases, we refrain from analyses based on their host prop- 
erties such as mass-weighted SFR or rest-frame color. 
One caveat of our luminosity-selected sample is a po- 
tential bias towa rds specific accretion modes such as the 
"Seyfert mode" ([Hopkins &i Her nquist 2006) or specific 
evolutionary stages in their fueling such as pre- or post- 
mergers. Fortu nately, we can test the later since merger- 
driven models (Hopkin s et al.ll2008al ) predict vastly dif- 
ferent properties (i.e., SFRs, morphology) of the hosts of 
AGNs before and after the coalescence of massive disk 
galaxies. 

We list in Table [1] the statistics of the parent galaxy 
sample and those hosting X-ray selected AGN. The red- 
shift and X-ray luminosity distribution of the full AGN 
sample (152) is shown in Figure Q] up to z = 1.02. Sym- 
bols denote their optical spectral properties with 82% 
lacking broad (FWHM > 1000 km s _1 ) emission lines 
and most characterized as narrow emission-line galax- 
ies. As further detailed in subsequent sections, we iso- 
late various subsamples based on redshift, galaxy mass, 
and AGN luminosity to optimize our analyses. For in- 
stance, as described above, we enforce an upper limit to 
the X-ray luminosity (log L% ax = 43.7) of the AGN when 
measuring a quantity with respect to its host galaxy. 
Also, we slightly increase our minimum X-ray luminos- 
ity (log L^ m =42.48) when determining the fraction of 
galaxies hosting AGN to ensure that a significant parent 
sample of zCOSMOS galaxies are capable of detecting 
each AGN to avoid effects based on limited statistics 21 . 

2.2. Stellar mass measurements 

Stellar masses, including rest-frame absolute magni- 
tudes (AB system; Mjj, My), are derived from fit- 
ting stel lar population syn t hesis models from the li- 
brary of iBruzual k. Charlotl (|2003h to both the broad- 
b and optical (CFHT : w, «, K s \ Subaru: B 7 V, g, r, i, 
z; ICapak et al.ll2007h and near -infrared (Spitzer /IRAC: 
3.6/x. 4.5/x: ISanders et al.ll2007h photometry using a chi- 
square minimization for each galaxy. The measurement 
of stellar mass (M») includes (1) the assumption of a 
Chabrier initial mass function, (2) a star formation his- 
tory with both a constant rate and an additional expo- 

21 Specifically, an overabundance of AGN at z ~ 0.35 falling 
at the flux limit of the JfMM-Newton observations, seen in Fig- 
ure [T] have an adverse effect on our determination of the fraction 
of galaxies hosting an AGN, due to the limited sample of zCOS- 
MOS galaxies capable of detecting AGN with Lx ~ 10 42 erg s — x , 
that disappears when implementing a slightly higher selection on 
luminosity. 
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Redshift 

Fig. 2. — Stellar-mass versus redshift for 7543 zCOSMOS galaxies 
(small grey circles). The box marks the mass-selected subsample of 
galaxies (log M > 10.6; 0.48 < z < 1.02) for which (specific) star 
formation rates are determined using [OIIJA3727. Galaxies hosting 
X-ray selected AGN are further marked by a larger black circle 
(filled: 42.0 < log L .5-10 kcV < 43.7; open: log L .5-10 kcV > 
43.7). 

nentially declining component covering a range of time 
scales (0. 1 < r < 30 Gy r ), (3) extinction (0 < A v < 3) 
following Cal zetti et al.l (|2000D . and (4) solar metallici- 
ti es. Further detail s on m ass measurements can be found 
in lBolzonella et all (|2008l ) and lMeneux et all (|2008h . For 
compatibility with the star formation rate calibration of 
iMoustakas et al.l (|2006[ ) implemented in subsequent sec- 
tions, we convert all masses to an equivalent based upon 
a Salpeter IMF by ap plying a multiplicative factor of 1.7 
(|Pozzetti et al J 120071) . We show the mass distribution 
for the parent galaxy population, including those harbor- 
ing X-ray selected AGN 22 with log L . 5-8.0 kcV > 42.0, 
as both a function of redshift (Figure [2]) and rest-frame 
color (U — V; Figure [3]). All masses are expressed in solar 
units (M Q ) throughout this work. 

We determine a minimum mass threshold that all 
galaxies must satisfy up to z ~ 1.0. The mass limit 
is set in order to ensure a fairly complete representation 
of both blue and red galaxies at all redshifts considered. 
In Figure El it is clearly evident that the mass limit of 
log M = 10.6 is essentially imposed by the red galaxy 
population at z > 0.8 (rightmost panel). The lack of red 
galaxies below this limit is due to our initial selection 
on apparent magnitude. iMeneux et al.l (|2008l ) estimate 
based on a series of mock catalogs from the Millennium 
simulation that the zCOSMOS 'Bright sample' is essen- 
tially complete for galaxies with log M w 10.6 at z = 0.8 

22 It is worth mentioning that there may be a potential problem 
that galaxies with even moderate-luminosity AGN may have inac- 
curate mass estimates; a bluer continuum will essentially reduce 
the stellar age and hence lower the mass measurements since the 
derived mass-to-light ratio depends strongly on the spectrum. 
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Fig. 3. — Rest-frame color U — V versus stellar mass split into three redshift intervals for galaxies and AGN with symbols described in 
Figure [2] The vertical line marks our imposed mass limit. 



while the completeness drops to ~ 50% at z = 1. In 
total, we have a sample of 2540 galaxies (0.1 < z < 1.02) 
above this mass limit of which 105 host X-ray selected 
AGN (see Tablc[T]for statistics regarding subsamples em- 
ployed herein.). 

3. [Oil] AS A STAR-FORMATION RATE 
INDICATOR AND THE AGN CONTRIBUTION 

Our primary aim is to use the emission line [OII]A3727 
to measure the SFRs of a well-defined sample of galaxies 
including those hosting AGN, as commonl y used for star- 
forming galaxies at these redshifts (e.g., iKewlev et al.l 
[200l iCooper et all I2001 . The availability of [Oil] 
within the spectral window of the zCOSMOS 'bright' 
program requires our sample of galaxies to fall within 
redshift range z = 0.48 - 1.02. The [Oil] line is one 
of a suite of spectral features (e.g., HaA6563, H/3A4861, 
[OIII]A5007) measured by an auto mated routine "plat- 
efitjvimos" (jLamareille et al.ll2008l ) that simultaneously 
fits all lines with a gaussian function to determine line 
flux and equivalent widths for the entire zCOSMOS 10k 
sample. Emission line fluxes are corrected for slit loss 
based on a comparison of their spectroscopic and pho- 
tometric (J-acs) magnitudes. Line measurements with 
a significance less than 1.15cr are quoted here as upper 
limits. 

An assessment of the AGN contribution to the [Oil] 
emission line is required even though the production of 
this lo w ionization l i ne has been shown to be relatively 
weak. ICroom et al.l (|2002D postulate, based on a com- 
parison of the equivalent width of [Oil] in the compos- 
ite 2dF QSO spectrum to that in a composite 'normal' 
galaxy spectrum, that [Oil] is mainly produced by star- 
formation in the host galaxies of AGN over a wide range 
of absolute magnitude. The [Oil] strength has been ob- 
serve d to be ~ 10 - 30% of the [O I II1A5007 emission line 
flux dFerland fc Osterbrockl Il986t iHol 120051 : iKim et all 
I2006H . iHol (|2005| ) has explored the feasibility of using 
[Oil] as an indicator of ongoing star formation for a sam- 
ple of PG quasars that in many cases provided only an 
upper limit. Upon fu rther investigatio n using emission 
line ratio diagnostics, IKim et ail (|2006f ) conclude based 
on a sample of ~ 3600 type 1 AGN (z < 0.3) selected 
from the SDSS that [Oil] emission is mainly attributed 
to AGN photoionization and does not have a strong HII 



component. 

The aforementioned results do not rule out the poten- 
tial effectiveness of using [Oil] as a SFR indicator for 
different samples of AGNs such those having higher lu- 
minosities, significant intrinsic obscuration or being at 
higher redshifts. For example, the stellar populations of 
'strong' type 2 AGN from the SDSS have young ag es sim- 
ilar to late-type galaxies (jKauffmann et al.ll2003tJ ) based 
on the age indicator £>„4000. We illustrate that this 
is likely related to the presence of ongoing star forma- 
tion. In Figure |4j we show the [OII]/[OIII] ratio for 
type 1 (panel a) and type 2 23 (panel 6;L[ ni] > 10 40 ' 58 
erg s- 1 ) AGNs from the SDSS having z < 0.3. The 
data for 25000 type 2 AGNs is taken from the high 
level products available from MPA based on the SDSS 
DR4 release; we require a S/N > 3 detection for [Oil], 
H/3, [OIII] and Ha. First of all, it is worth highlight- 
ing that the best-fit linear relation to the type 1 AGN 
(panel a) illustrates a scenario where [Oil] emission scales 
with [OIII] (log i[on] tx 0.36 x £[0111]) but not with 
a one-to-one correspondence. Even without correcting 
for extinction, we see that there is evidence (Fig. \4$>) 
for elevated [Oil] /[OIII] for type 2 AGNs as seen by the 
flatter slope of the best-fit relation (slanted dashed line; 
log -£z[on] oc 0.58 x Lronil) compared to that of the type 
1 AGNs. The change in slope appears to signify that the 
strength of an additional component to [Oil] (i.e., star 
formation) increases with AGN luminos i ty po ssibly re- 
lated to the results of iKauffmann et"afl (|2007l ) where a 
decline in stellar age is seen for AGNs with higher mass 
accretion rates. Finally, we see conclusively that the 
AGN host galaxies in zCOSMOS, while having higher 
[OIII] luminosities, have elevated [Oil] /[OIII] ratios com- 
pared to type 1 AGNs (panel a), and exhibit a similar 
slope to the type 2 AGNs in the SDSS and slightly fur- 
ther enhanced [Oil] emission. 

When correcting for internal extinction based on the 
Balmer decrements, we find that the type 2 AGNs have a 

23 The luminosity selection of type 2 AGNs from the SDSS is 
eq uivalent to the de f inition of 'strong' AGN (L > 10 7 Lq) given 
in IKauffmann et alj d200 3b). This luminosity selection, based on 
extinction-corrected values, guarantees that the sample is domi- 
nated by Seyfert galaxies comparable to the luminosities of the 
zCOSMOS AGN. 
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Fig. 4. — Emission-line properties of AGN. (a) Observed (i.e., no extinction correction) line ratio [OII]/[OIII] versus [OIII] luminosity for 
zCOSMOS galaxies (log M > 10.6; small green poi nts) and those ho sting AGN shown by larger red circles. For comparison, emission line 
properties of type 1 AGN from the SDSS (z < 0.3; |K im et al. 2006) are marked by the small black dots that are further characterized by 
the best fit linear relation (solid line) and mean ratio < [OII/[OIII] >= 0.27 (open white square). The dashed horizontal line denotes our 
single assumption for the value of a purely A GN dominated [QUI/ [011 1] ratio (see text for further details), (b) Same as panel a but with 
type 1 AGNs replaced by type 2 AGNs from [Kauffmarm et al. (2003b). The slanted lines are the best fit linear relation to type 1 (solid) 
and type 2 (dashed) AGNs. (c) Histogram of the extinction-corrected [OII]/[OIII] distribution for the type 1/2 AGN (short dash/long 
dash) from the SDSS and our zCOSMOS sample (solid) all having log £[oiii] > 41.3. 



significantly higher distribution of [Oil] /[OIII] than the 
type Is (Fig. 0J;). This further suggests that star for- 
mation is more prevalant in the type 2 AGN popula- 
tion compared to those with an observable br oad line re- 
gion; s uch a hypothesis has been put forth bv lKim et al.l 
(2006) to explain the [Oil] emission (< [OII]/[OIII] >= 
—0.12) evident i n the more luminous, type 2 QSOs 
(|Zakamska et al.ll2003D . To further justify our use of 
[Oil] , we demonstrate in Figure [5] that a correlation be- 
tween total [Oil] emission (HII+AGN) and the spectral 
index £>„4000 exists and is in agreement with significant 
levels of ongoing star formation seen in t he hosts of type 
2 AGN (|Gu et al.ll2l)0irYan et al.ll200l . 

For our purpose, we use the higher ionization line 
[OIII]A5007 when present in our spectra, and an empir- 
ical relation between [Oil] and [OIII] to statistically re- 
move the component of the detected [Oil] line that can be 
attributed to the AGN. This entails a n assumption that 
the [Q UI] line is purely of AGN origin (flCauffman n et all 
2003b) that may lead us to underestimate SFRs due to 
a significant stellar contribution but only in very few 
cases since the [OIII] luminosities of Seyfert galaxies in 
SDSS, and also zCOSMOS AGN hosts (see Figure gb), 
are substantially highe r than that of typi cal HII galax- 
ies (|Kewlev et al.ll2006l) . iKim et al.l (2006) find that the 
median value of the observed [Oil] /[OIII] ratio is 0.27 
with significant dispersion (see Fig. 0Ji). We chose to 
use a slightly lower ratio ([OII]/[OIII]=0.21) that is the 
mean value for type 1 AGN in the SDSS sample (M. 
Kim private communication) with log Liqiii] > 41.5, a 
luminosity regime similar to our zCOSMOS AGN sam- 
ple. This higher luminosity cut guarantees that the ratio 
best reflects that produced by the AGN, while effectively 
minimizing a contribution from HII regions. This results 
in slightly higher SFRs for our AGN sample than if we 
chose to use the mean of the entire type 1 SDSS sample 
but our final results in this study are consistent using 



either value. We make no attempt to use a luminosity- 
dependent correction, as inferred by the SDSS type 1 
AGN, since a simple linear relation may not be evident 
at higher luminosities (log L\oiii\ > 42). Also, there is 
no need to correct this relation for extinctio n given the 
low le vels of dust attenuation (Ay « 0.2 mag: IKim et al.l 
l2006f ) observed in these type 1 AGN. We strongly note, 
given the high dispersion in the [Oil] /[OIII] ratio for type 
1 SDSS AGN, that this is only a statistical correction ap- 
plied to the population as a whole to infer the global or 
mean properties of the sample; the individual measure- 
ments of [Oil] strength associated with star formation in 
a particular galaxy are expected to be inaccurate. We 
further note that even with this correction the more lu- 
minous AGN tend to have higher SFRs as evident by the 
best-fit relation for type 2 AGNs in the SDSS (Fig [4f>; 
log ipii] c< 0.37 x ijoiii])- I n addition, we assess the 
impact of the dispersion in the [Oil] /[OIII] ratio for type 
1 SDSS AGNs by assuming a normal distribution of the 
line ratio (< log[OII]/[OIII] >= -0.69; a = 0.25) and 
perform many iterations in our determination the distri- 
bution of SFR. 

We consider extinction due to dust as an important fac- 
tor in determining the component of the observed [Oil] 
emission attributed to AGN photoionization given that 
our sample is predominantly composed of type 2 AGN 
(see Fig. [TJ. An extinction-corrected [OIII] luminosity 
is found for each zCOSMOS AGN with detected [OIII] 
line emission using Ay — 0.8. Since most of our sample 
does not have Ha and H/3 within our observed spectral 
window, we chose to implement a level of extinction 24 
based on the mean Balmer de crement of type 2 Seyferts 
in SDSS (|Kewlev et al.ll2006l ). This same level of at- 

24 This amount of optical attenuation is similar to that of X-ray 
selected AGN; we find a median Ay of 0.98 based on the Balmer 
decrements of 20 AGN in our sample with significant Ha and H/3 
line measurements (Mainieri et al. in preparation). 



Star formation in AGN host galaxies 



7 



42.5 
42.0 

T 41.5- 
41 .0 

o 

I 

cn 40.5 - 

40.0 

39.5 
1 .0 





. 1 1 '■' 




ii 

(a); 


- - 










■ ' 






- ■■ 
























;-'.^;.'V.V.v>'.'>y.^ 






i 







' .2 



1.4 1.6 

D„4000 



2.0 



.4;§£^2.^i : ' 



cn 
O 



0.5 



o.o - 



-0.5 - 



■'" ■}:-} .*.V.</i'Sf*::?'.;--'. : 



1 .2 



.4 1 .6 
D„4000 



1.8 



2.0 



Fig. 5. — Relation between D n 4000 and (a) [Oil] luminosity and 
(b) [OII]/[OIII] ratio for 4609 luminous(L [om > 10 7 Lq), type 2 
AGN from the SDSS. In panel ft, the median [Oil] /[OIII] ratio is 
shown (large red circles) in bins of D n 4000 having a width of 0.1. 

tenuation is then reapplied to the inferred AGN compo- 
nent to [Oil], based on the empirical relation mentioned 
above, and subtracted from the observed (i.e. no dust 
correction) [Oil] emission line luminosity to provide an 
estimate free of any AGN contribution. We are confi- 
dent that this method is applicable to the zCOSMOS 
sample given that the [Oil] luminosities of AGN hosts 
are all systematically higher (~ 3 — 4x) than that ex- 
pected from gas photoionized by an AGN as indicated by 
the SDSS type 1 AGN (Fig. \%jp) of similar luminosities 
(log i[om] > 41.3). Also, we highlight that the enhance- 
ment of [Oil] relative to [OIII] for our zCOSMOS AGN 
is not purely induced by our extinction corrections as ev- 
ident in observed relation (i.e., no extinction correction; 

Fig. EM- 

For AGN host galaxies at z > 0.8, we need to modify 
our method since [OIII] is no longer within our observed 
spectral window. Fortunately, we can utilize the strong 
correlatio n between hard X-ray ( 2-10 keV) and [ QUI] lu- 
minosity (H eckman et al.l l2005t iPanessa et al.l I2006D to 
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Fig. 6. — (a) Extinction-corrected [OIIIJA5007 luminosity versus 
X-ray luminosity (2-10 keV) for AGN. A clear correlation is evident 
as shown by the best-fit linear relation to these da ta (solid line) 
and local Seyferts (dotted line; Panessa et al. 2006). (ft) Fraction 
of the [Oil] line luminosity attributed to an AGN and split by the 
method of determining the AGN contribution. The number distri- 
bution of galaxies with detected [OIII] emission, used to estimate 
the amount of [Oil] due to AGN activity, is given by the solid his- 
togram. Galaxies without [OIII] detections are indirectly assessed 
by their hard X-ray luminosity (dashed histogram). There is no 
evidence for a systematic offset between both methods. 



estimate the AGN contribution for these cases. It is im- 
portant to highlight that our AGN are X-ray selected by 
the 0.5-10 keV band thus limiting the inclusion of more 
heavily abs orbed type 2 AGN such as the Compton-thick 
population (|Fiore et al.ll2008h . for which this correlation 
is less evident if the X-ray luminosities are not corrected 
for X-ray absorption (|Heckman et aLlfeOOSD as is the case 
here. In Figure [5Jj, we plot the [OIII] line luminosity, 
corrected for extinction as done in the previous para- 
graph, as a function of the rest-frame 2-10 keV luminos- 
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ity for AGN in the zCOSMOS sample with redshifts over 
a larger baseline (0.2 < z < 1.0). We use these lower red- 
shift AGN in order to provide higher statistics to measure 
the best-fit linear relation between these two quantities. 
We determine the best fit Lx — L^jiii] relation (Equa- 
tion [T]) for our sample by implementing a bivariate linear 
regression (EM) algorithm with ASURV (Survival Anal- 
ysis fo r Astronomy package Rev. 1.2; lLaVallev et al.l 
|1992() that considers limits as well as detections. 



logL [OIII] = (0.729±0.101) log L 2 -io fce v + (10.307±4.381) 

(1) 

This relation (Fig. [6ji), as shown by the solid line, is 
similar to previ ous studies based on a sample of low red- 
shift AGN fe.g. JPanessa et al.l2006f ) shown by the dotted 
line. We find a nearly equivalent slope although a higher 
normalization (i.e., larger values of £[oni] for & given 
£2-10 kcv)- The strong linear relation (Pearson corre- 
lation coefficient r = 0.63) strengthens our assumption 
that most of the [OIII] line luminosity is due to AGN 
photoionization. Any significant contribution from stars 
to the [OIII] luminosity would increase our derived star 
formation rates in AGN hosts. In Figure O), we demon- 
strate that there are no systematic offsets by introduc- 
ing this indirect probe of the [OIII] line luminosity; both 
methods essentially remove a similar fraction (< 60%) of 
the [Oil] luminosity attributed to AGN photoionization. 

We determine S FRs u sing the empirical calibration of 
iMoustakas et al.l |2006). This relation considers inter- 
dependent factors (i.e., dust extinction, metallicity, and 
ionization) to provide a SFR tracer that is applicable 
to samples for which these quantities are inadequately 
known. To derive the following empirical relation, we fit 
a linear relation to the bright end (log L(B) > 9.5) of 
the dependence of SF R on absolute .B-band magnitude 
(M B ; see Figure 19 of IMoustakas et al.ll2006| ): 

log SFR{[Oll}) =log i [O n]-41-0.195xM B -3.434 (2) 

The AGN component to the emission line luminosity 
L\oii] is removed as described above and no correction 
for dust extinction is applied that is in essence an as- 
sumption that star formation is external to a dusty NLR 
and circumvents the potential problem that extinction 
of t he nuclear regio n may differ from that of HII regions 
(see lGu et al.l 12006). We do not attempt to adjust (i.e., 
lower) the upper limits on nondetections by considering 
the contribution of an AGN. It is worth noting that an 
additional caveat of our method is that there is an un- 
derlying assumption that the metallicity and extinction 
of galaxies with or without AGN are similar. 

4. PHYSICAL PROPERTIES OF AGN HOST 
GALAXIES 

4.1. Stellar masses 

We measure the fraction of galaxies that host X-ray 
selected AGN as a function of their stellar mass to deter- 
mine the characteristic mass above which AGN activity is 
most prevalent. First of all, it is apparent in FigureOthat 
galaxies hosting X-ray selected AGN are preferentially 
massive (log M > 10.5) as evident over the full redshift 
range and in agreement with most studies to date (e.g., 




Fig. 7. — Stellar mass distribution of galaxies hosting AGN in 
the redshift range 0.1 < z < 0.5. Fifty-two galaxies {log M* > 
9.5) with X-ray selected AGN (42 < Log Lo. 5 _io keV < 43.7) are 
shown by the solid histogram. The observed parent distribution 
of 3356 galaxies (dotted histogram) is renormalized to match the 
AGN distribution. The data points represent the AGN fraction 
(42.48 < Log -£<o.5-iO kcV < 43.7; scale on the right hand vertical 
axis) in fixed mass intervals (shown by the horizontal bars centered 
on each data point) with la error bars. 



iKauffmann et al.ll2003bHBundv et al.ll2008f ). For this ex- 
ercise, we consider the redshift interval 0.1 < z < 0.5 for 
which we have a high degree of completeness for galax- 
ies over a wide mass baseline (log M > 9.5). We follow 
the technique discussed in § 3.1 of iLehmer et alJ (2007) 
to determine the AGN fraction of our parent popula- 
tion of galaxies that accounts for the spatially varying 
sensitivity limits of the X MM observations of the COS- 
MOS field (see Fig. 17 of iCappelluti et al.l [2007h . The 
necessity of this approach is demonstrated in Figure [TJ 
which shows the limiting X-ray luminosity as a function 
of redshift for the entire galaxy sample. Even though 
the sensitivity of the XMM coverage is remarkably uni- 
form, some dispersion is present as shown by the rela- 
tively narrow distribution of the X-ray upper limits at 
each redshift. To properly account for the luminosity- 
rcdshift relation, we determine the contribution of each 
AGN separately to the total fraction. The AGN fraction 
(/; see equation[3]below) is determined by summing over 
the full sample of AGN (N) with Agai^ representing the 
number of galaxies in which we could have detected an 
AGN with X-ray luminosity L^- The different sampling 
rates, based soley on slit placement, of the galaxies (S ga i) 
and AGN (i.e., X-ray sources 25 , S x ) are taken into ac- 
count in equation[3l We estimate the associated la error 
(equation |4|) using binomial statistics where N^ g l is the 
number of AGN that would be detected if all galaxies 
have the same limiting X-ray sensitivity and the sample 
of AGN was randomly selected. Here, we only consider 
AGN with 42.48 < log Lo.s-iokeV < 43.7. As previously 
mentioned, the lower limit ensures that we have a star 

25 In many cases, X-ray sources have redshifts based on a random 
placement of a slit for which the sampling reflects the overall galaxy 
population (Sx = S ga i) 
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tistically significant sample of parent galaxies (> 700) 
that could host each AGN while the upper limit restricts 
the sample to low-to-moderate luminosities thus securing 
the accu racy of their hos t -galaxy masses. We refer the 
reader to lSilverman et a l. (2008a'|) for further details and 
results employing this method. 

t = V ^ (3) 

1 t[ N salti /S gal { } 

^ = K{1 >< (Ngal ~ Kff)/N 3 gal (4) 

In Figure we plot the fraction of galaxies that host 
AGN and the number distribution (TableU Sample A) of 
both galaxies and those with AGN over a slightly wider 
luminosity interval (42.0 < log £o.5-iokcV < 43.7). The 
absolute fraction is lo w (~ 1 — 4%) due to our re strictions 
on X-ray luminosity. Kauffmann et al. (2003b) measure 
a fraction ~ 5x higher although they include in their 
sample AGN ~ 80 x fainter. The difference in these 
fractions is mainly due to the steep fa int-end slope of 
the lu minosity function at low redshifts (fShinozak i et alj 
2006). Here, we are mainly interested in the relative 
change of the AGN fraction with galaxy mass. Clearly, 
we see that the fraction of galaxies that host AGN mono- 
tonically increases with stellar mass simila r to results 
based on obscured dKauffmann et aT]|2003bO and radio- 
loud (|Best et al.l l2005t) AGN in the SPSS. 

4.2. Star formation rates 

In Figure H we show the SFRs for 1820 zCOS- 
MOS galaxies (Table [Q Sample B 26 ) including 43 
of them that host moderate- luminosity AGN (42.0 < 
log Lo. 5-10.0 koV < 43.7). The results are presented be- 
fore (panels a and b) and after (panels c and d) the 
correction for the AGN contribution to the [Oil] lumi- 
nosity based on a single-value (0.21) of the observed, 
median [OII]/[OIII] ratio of type 1 AGNs in the SDSS. 
First of all, we see that 88% of AGNs have detectable 
[Oil] emission 2 (panel a). It is apparent by comparing 
the SFR distribution for AGN hosts in panels b and d 
that a significant shift in the distribution occurs when 
removing the AGN contribution to the [Oil] emission 
line. Even so, we find that the SFRs of AGN host 
galaxies are almost exclusively between 1-100 Mq yr" 1 
(Fig. [5fc,<i), a range consistent with analytic models o f 
AGN hosts with star-forming disks (|Ballantvnel l2008f) 
and sup portive of constraints on the cosmic IR back- 
ground (jBallantvne fc Papovichll2007t) . 

To quantitatively test whether the SFR distribution 
of AGN hosts differs from that of the underlying galaxy 
population, we have determined the SFR distribution, 
including those with upper limits, while also considering 
the dispersion in the AGN contribution to the observed 
[Oil] emission-line luminosity. We assume a log nor- 
mal distribution of the [OII]/[OIII] ratio with parameters 

26 The sample is slightly smaller than that given in Table [T] four 
AGNs have [Oil] falling outside the observed spectral window that 
can vary slightly from slit-to-slit. 

27 The fraction of AGN with detected [Oil] emission is 68% 
(S/N > 2) and 61% [S/N > 3) depending on the given line signifi- 
cance. We have confirmed our results based on the lower S/N ratio 
(1.15) with those based on these smaller, high confidence samples 
of AGN. 



given in Section [3j Based on one thousand iterations, we 
determine the mean SFR distribution and perform sta- 
tistical tests on each individual iteration compared to the 
parent galaxy population. The results of this exercise are 
shown in Figure [5] for two AGN samples of differing lu- 
minosities. First, we consider AGN, as done above, that 
fall within the luminosity interval 42 < log Lx < 43.7 
for which we have high confidence in our well-matched, 
mass-selected parent sample of galaxies. Based on this 
limited sample of AGNs, we cannot significantly discrim- 
inate between the two distributions (Fig. \9]a) given the 
results of the individual K-S tests (panel b). Although, 
we do find that the median SFRs shown in panel c are for 
the most part higher than that of the full galaxy sample 
(4.7 Mq yr -1 ) in each iteration. 

We are able to use a less restrictive sample that in- 
cludes higher luminosity AGN (log Lx > 43.7) thus im- 
proving our statistics although at the expense of main- 
taining reliable stellar masses across our sample. In Fig- 
ure [SJi — /, we present the equivalent analysis for AGNs 
with log Lx > 42. Based on this larger sample, we 
see that the SFR distribution of AGN hosts is shifted 
to higher values than the overall population (panel d). 
The K-S tests on the individual distributions now re- 
ject the null hypothesis (i.e., distributions are equiva- 
lent) at the ~ 99% level (> 2.5a) in essentially all it- 
erations (panel e). Moreover, the median SFRs are all 
systematically higher (SFR ~ 9.5 M Q yr -1 ) than that 
of the full galaxy sample. We further point out that 
the SFR distribution of AGN hosts while being elevated 
in comparison to the underlying massive galaxy popula- 
tion is essentially equivalent to those forming stars (i.e., 
emission- line galaxies). Therefore, we conclude that a 
plentiful gas reservoir is a necessary ingredient for the 
fueling of AGN as indicated by the presence of signifi- 
cant star formation with rates reaching up to ~ 100 M Q 
yr" 1 . 

We investigate how star formation in galaxies hosting 
AGN is evolving compared to that of the parent galaxy 
population. In Figure [TU1 we show SFR as a function of 
redshift for galaxies with log M > 10.6 and those host- 
ing AGN (Table [2 Sample B). A systematic shift of the 
distribution of zCOSMOS galaxies towards higher SFRs 
with increasing redshift is clearly evident with those host- 
ing AGN exhibiting a similar behavior. Therefore, we 
find that the SFRs of AGN hosts are dictated by that of 
the underlying galaxy population thus solidifying simi- 
lar evidence base d on the color evolution of AGN hosts 
in the E-CDF-S (|Silverman et alj l2008al) . As a conse- 
quence, ongoing star formation is most likely responsi- 
ble for the significant population of AGN hosts over the 
reds hift range 0.5 < z < 1-4 having blue rest-frame col- 
ors dSanchez et al.lfeooj iBohm et al]|2007t iNandra et all 
l200l " with evidence for such a trend remaining in place 
for quasar hosts at higher redshifts 1.8 < z < 2.8 
(jJahnke et al.ll2004bt) . 

It is illuminating to compare the SFRs from zCOS- 
MOS to AGN at lower redshifts. For consistency 
with zCOSMOS, we me asure SFRs of the type 2 AGN 
from SDSS (z < 0.3; iKauffmann eTall l2003bl) using 
the [ OII]A3727 and [O I II]A50 07 emission line fluxes 
from [Brinchm ann et all (|2004f ) to determine an AGN- 
corrected [Oil] luminosity. We further select only those 
galaxies having stellar masses above our zCOSMOS limit 
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(d) Minus AGN contribution 
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Fig. 8. — Star formation rates (Mq yr _1 ) versus stellar mass for galaxies (0.48 < z < 1.02; log Al* > 10.6) before (a) and after 
(c) the removal of the AGN contribution. Measurements are shown by either a small black (emission-line galaxy) or large red (AGN; 
42.0 < log -Lrj.5— 10 keV < 43.7) circle. Upper limits are shown as an arrow with their color descriptive of the type as given above, (b, d) For 
ease of comparison, the distributions, including those with upper limits, are shown by either a solid (AGN hosts) or dashed (all galaxies) 
histogram with the full galaxy sample scaled down to match the AGNs. 



{log M* > 10.6) based on the conversion factor ( Kroup a 
to Salpeter IMF) given in iBrinchmann et al I d2 004) . 
Here, we use the calibration of iKewlev et al.l (|2004l ) to 
derive SFRs 28 in order to have consistenc y with the 
SFRs measured for the type 1 SDSS sample (|Kim et al.1 
2006). Initially, we demonstrate that the ongoing SFRs 
sh own in F igure 1101 are in a greement with the findings 
of iKauffmann et al J (l2003bh based on L>„(4000): weak 
AGNs (i[o//7] < 10 '1/0 ; blue points) have low SFRs 
(-0.2 M Q yr- 1 ) while strong AGN (i [0 777] > 10 7 L Q ; 
small red points) are more actively forming stars (~1 
M Q yr -1 ). For ease of visualization, we show the best- 
fit relation log SFR oc log(l + z) for only the zCOS- 
MOS galaxies with significant [Oil] emission (upper lim- 
its are not considered; black curve) and those hosting 
AGN (red curve). Remarkably, we further find that the 
SDSS AGN may be the low redshift analogs of the AGNs 

28 We have confirmed that the Moustakas et al. (2006) relation 
gives the same results. 



in the zCOSMOS survey given their close proximity to 
an extrapolation of the evolution of zCOSMOS galaxies. 
In addition, the low-to-moderate levels of star formation 
(« 0.5-3 M yr" 1 ) in type 1 SPSS AGN (shown rou ghly 
by the green triangle in Fig. [lOj iKim etall I2006T) are 
in agreement with the type 2 AGN from SDSS and the 
aforementioned passive evolutionary scenario. In light of 
these results from the SDSS that effectively extend our 
redshift baseline, we reiterate our conclusion that the 
SFRs of AGN host galaxies are reflective of the overall 
star formation history of galaxies and provide no indica- 
tion of the suppression or truncation due to a mechanism 
related to the AGN itself. 

4.3. Stellar ages 

To complement our study of the ongoing SFRs of AGN 
hosts , we use the spectral index D„(4000) (jBalogh et al.l 
Il999h . determined by our " platefit_vimos" routine for 
each galaxy in the zCOSMOS sample to infer the age of 
the overall stellar population on longer timescales (> 0.1 
Gyr). This index is the ratio of the average flux density 
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Fig. 9. — Star formation rates of AGN hosts with dispersion in the applied correction for the AGN contribution. The two rows are given 
for AGN spanning a different luminosity range: (a — c) 42 < log Lx < 43.7, (d — f) log Lx > 42. The left panels (a, d), equivalent to 
those in Figure [5] show the mean distribution of SFR, including upper limits, based on 1000 iterations. The middle panels (ft, e) give the 
probability distribution based on a KS tests that the two samples are equivalent. The right panels (c, d) display the distribution of the 
median SFR for each iteration with the median SFR (4.73 Mq yr _1 ) of all galaxies depicted by the vertical dashed line. 



F x in the continuum bands 3850-3950 A and 4000-4100 
A. This is essentially a measure of the strength of the 
4000-A break with galaxies having experienced a recent 
episode of star formation exhibiting a smaller index due 
to the presence of young stars. In Figure [TTl we show the 
values of D„(4000) for galaxies with respect to their spe- 
cific SFR. We implement two mass limits (log M > 10.6: 
panels a, b; log M > 11.1: panels c, d) in order to check 
for consistency when the underlying £>„(4000) distribu- 
tion of galaxies is significantly different: the higher mass 
cut results in a distribution dominated by more evolved 
galaxies as evident by its peak at Z3„(4000) ~ 1.8. Ini- 
tially, it is apparent (Fig. [Tib ) that there is a good cor- 
respondence between Z3„(4000) and specific SFR for all 
galaxies includi ng those with AGNs. We then use the 
relation given in iKauffman n et all (|2003aj ) that assumes 
an instantaneous burst model of solar-metallicity to infer 
stellar age from the value of D at (4000) as given by the 
scale bar in Figure [TTa: this relation provides us with a 
rough assessment of the actual ages since many of the 
zCOSMOS galaxies including those with AGN may have 
had a more sedative existence in their recent past. As a 
result, we find that most galaxies with AGN contain a 
young stellar component since 70% have D„(4000) < 1.6 
( a 9 e ^5 2 Gyr; Fig Ilia , c). Furthermore, we measure 
the fraction of galaxies hosting an AGN as a function of 
£>„(4000) and demonstrate (Fig. [Tib, d) that a galaxy is 



more likely to have an accreting SMBH if there is a suf- 
ficient supply of gas, fully consis tent with the mass mea- 
surem ents of HI (|Ho et al.ll2008[ ) and CO (|Scoville et all 
l2003f ) in AGN hosts, given the higher rate of occurance 
in galaxies with younger ages. These results effectively 
extend such studies base d on D„(4000) at low redshift 
(jKauffmann et aLll2003bl ) up to z ~ 1. 

There is a noteable discrepancy between our claim 
for a higher level of AGN activity for star-forming 
galaxies and the location of AGN hosts on the color- 
magnitude diagram. Both X -ray selected (jNandra et al.l 
l2007tlSilverman et al.l2008aft and optically-selected AGN 
(|Martin et al.l [20071) exhibit a low AGN fraction along 
the red sequence, enhanced activity in the intermedi- 
ate region (i.e., "green valley") between the blue and 
red galaxy populations, and a dropoff towards very blue 
galax ies. On the contrary, the mass-selected AGN sam- 
ple of lKauffmann et aTl (|2003bh clearly have stellar prop- 
erties very similar to late-type galaxies (see their Fig. 14) 
and, thus a low fraction of AGN within the 'blue cloud' 
reported b y the aforemen t ioned studies is surprising. In- 
terestingly, [Lehmer^etlal] ((2008) find that the fraction of 
the stacked X-ray signal of late-type galaxies attributed 
to AGN emission continuously rises with SFR. 

To investigate this further, we measure the fraction of 
galaxies hosting AGN as a function of their rest-frame 
optical color (U — V) for both a luminosity and mass- 
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Fig. 10. — Cosmic evolution of star formation. At z > 0.48, we show the SFR-2 distribution for all zCOSMOS galaxies with log M > 10.6 
(small black circles and crosses) and those hosting AGN with log Lx > 42 (73; large red circles and arrows). The best-fit linear relation 
for the emission-line objects is shown for both zCOSMOS populations (black: galaxies; red: AGN hosts) with an extrapolation to lower 
redshifts (dashed lines). For compa rison, we plot SFRs of AG N hosts from the SDSS with an equivalent selection on stellar mass; obscured 
AGN (type 2) from the sample of iKauff mann et al. (2003Q) are shown with strong AGN (log Lq jjj > 40.5) in re d and those of lower 
luminosity in blue. A large green triangle marks the mean value of the SFR for SDSS type 1 AGN l)Kim et~aUl2006T) . 



selected sample (Figure fT2h. b). Since we are not con- 
strained by the use of [Oil] for this exercise, we have 
extended the redshift baseline 0.1 < z < 1.02 to im- 
prove the statistics. In the first panel, we find that 
the fraction of galaxies with X-ray selected AGN is 
strongly peaked in the " green valley" for the luminosity- 
selected sample ( panel a ) thus agreeing with the results of 
iSilverman et al.l ((2008a). In contrast, the mass-selected 
sample, behaves similarly for colors U — V > 1.4, but 
does not have a decline towards bluer colors U— V < 1.4. 
This is easily understood since blue galaxies are known 
to have lower mass-to-light ratios and the AGN frac- 
tion rises substantially with host galaxy mass as demon- 
strated in Section [4~T1 We conclude that blue, star form- 
ing galaxies do have enhanced levels of AGN activity 
and there is no evidence for a significant delay in the 
emergence of nuclear activity with respect to the onset 
of star formati on. This result is consis tent with the re- 
cent findings of lSilverman et aT] <|2008af l . based the mor- 
phology of the host galaxy, that blue (U — V < 0.7), 
bulge-dominated (n sers i C > 2.5) galaxies have the high- 
est incidence (21.3%) of AGN activity, and the blue rest- 
frame colors of low redshif t (z < 0.2) quasars irrespec- 
tive of their morphology (|Jahnke et al.l l2004af ) . How- 



ever, there remains the possibility that the lower frac- 
tion of red galaxies hosting AGN is due to some form of 
self-regulating feedback that effectively reduces the AGN 
luminosity below our flux limit. Nonetheless, these re- 
sults lend no support for a simple model prescription 
that attributes the truncation of star formation or the 
color evolution of galaxies to AGN feedback and bring 
into question whether studies based on the optical emis- 
sion line diagostics of AGN activity miss a significant 
number of those resid ing in star-forming galaxies (see 
Schawins ki et al.lfeOOTl for a discussion on this topic). 

5. AGN-STAR FORMATION CONNECTION AND 
CO-EVOLUTION 

Having established an association between AGN ac- 
tivity and star formation, we are motivated to deter- 
mine how closely these phenomena are related on a case- 
by-case basis that can potentially signify an underly- 
ing causal connection. The existence of such a rela- 
tionship may be realized given the recent findings that 
AGN accretion power as probed by [OIII] luminosity 
is higher for galaxies w ith younger stellar populations 
(|Kauffmann et aLll2007f ). Ideally, we would like to inves- 
tigate such a relation using a quantity more closely as- 
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Fig. 11. — (a) Specific SFR versus the stellar age ind i cator D n (4000). Symbols are equivalent to those in Figure |8ji. An approximate 
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sociated with the accretion process, namely X-ray emis- 
sion, known to originate closer to the black hole and an 
indicator of the instantaneous star formation rate. 

To do so, we have plot in Figure [T3k the SFRs of our 
AGN sample, as measured by [Oil] strength, versus hard 
X-ray luminosity (2-10 keV) and inferred mass accretion 
ra te while implementing a bolometric correction as given 
in iMarconi et al.l (|2004l ) and an accretion efficiency of 
0.1. We find that a weak correlation exists based upon a 
Pearson correlation coefficient of 0.17 and a linear fit that 
has a shallow slope with significant dispersion in its value 
(0.28 ± 0.22). We conclude that underlying complexities 
such as the efficiency of transferring gas to nuclear re- 
gion over kiloparsec scales, and varying duty cycles for 
star formation and accretion may contribute to the large 
dispersion in these relations. 

There have been some recent claims that an obscured 
phase, coupled with enhanced star formation, may rep- 
resen t an early stage in the subsequent evolut i on of AGN 
(e.g. JPage et all 12004 lAlexander et al.ll2005t iKim et all 



120061 iHopkins et aT1l2008at lPollettall2008h . To test this 
scenario, we have marked those AGN in Figure fT3k that 
have excessive X-ray absorption (Nh > 10 22 cm -2 ) based 
on their hardness ratio [HR ={H - S)/H + S)> -0.2] 
determined by the X-ray counts in the soft (S: 0.5-2.0 
keV) and hard (H: 2-10 keV) bands assuming a pow- 
erlaw spectrum with a photon index of 1.9 at an ef- 
fective redshift of z ~ 0.7. We find that the absorbed 
sources span the same range of SFR as the unabsorbed 
sources. A KS test give a probability of 65% that both 
absorbed and unabsorbed AGN could be drawn from the 
same parent population. These results are not depen- 
dent on the chosen division in hardness ratio. Therefore, 
we conclude that an amendment, such as the aforemen- 
tioned evolution ary scenario, to the unification model 
(|Antonuccilll993l ) for these moderate-luminosity AGNs is 
not supported by our findings. Although, there remains 
substantive evidence that enhanced star formation may 
be associated with n uclear obscuration for the more lu- 
minous QSOs (e.g.. iPage et al.ll2004 lLacv et alj 120071 : 
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Fig. 12. — Rest-frame color distribution U — V of galaxies (0.1 < z < 1.02) and the fraction hosting AGN. We show the distribution for 
58 galaxies (solid histogram) selected by optical luminosity (My < —21.5; panel a) and 65 galaxies above a fixed mass limit (log M > 10.6; 
panel b) that host AGN. For comparison, the distribution for the parent galaxy population (dashed histogram) is normalized to match the 
AGNs in each panel. The decline in the AGN fraction towards the bluest colors (panel a) appears to be due to the inclusion of galaxies 
with low mass-to-light ratios that are not present in the mass-selected sample (panel b). 



IZakamska et al.ll2008f) possibly due to different physical 
mechanisms for triggering mass accretion on SMBHs. 

It is of much interest to determine the relative growth 
rate of a galaxy and its central SMBH as a function 
of redshift in light of the well- established local relation 
MzK/Mhdaz ~ !-5x 10" 3 (e.g- lMcLure fc Dunlorj|200l 
lHaering fc Rixll2004D . We plot in Figure [Hb the ratio 
of mass accretion rate onto a SMBH to the SFR as a 
function of redshift. The median value (1.9 x 10~ 2 ) is 
roughly an order-of-magnitude higher than the local ratio 
Mbh /Mbuige ■ This difference is likely due to the varying 
timescales between SMBH accretion and star formation 
M BH /At gal a Mbh /Mbuige X SFR/dM accr dt^ a 0.1. 
We note that our SFR measurements most likely in- 
clude a significant disk component thus a more rigor- 
ous assessment is required and beyond the scope of this 
work. Nonetheless, by assuming that star formation oc- 
curs over rougly a dynamical timescale At ga i a 10 9 yr, 
the duty cycle for AGN activity (a 10 8 yr) is consis- 
tent with the luminosi ty-dependent model predictions of 
iHopkins et all (|2008bf) for SMBHs with M BH ~ lfj 8 M 
and accreting above an Eddington ratio of 0.01, both 
within a physical regime spanned by our sample assum- 
ing these AGNs have already settled on a SMBH - bulge 
relation at their respective redshifts. Alternatively, if 
we assume a timescale for galaxy growth to be the in- 
verse of their sSFRs as shown in Figure QTJ the AGN 
lifetimes can reach up to ~ 10 10 yr that fully illustrates 
that these AGNs and their host galaxies are close to be- 
ing fully matured although have SMBH growth times 
at leas t an order-of-magnit ude less than those in the 
SDSS (|Heckman et al. I [2004( 1. Furthermore, we find that 
there is no dependence of the ratio {dM accr dt~ l / SFR) 
on redshift thus supporting a co-evolutionary scenario 
where both the average SFR and mass accretion rates 
onto SMBHs are rapidly declining with equivalent rates 
from z <~ 1 to the present possibly due to diminishing 
fuel supplies. 



These comparisons suggest that beyond the local uni- 
verse the Mbh — Mbuige relation should display higher 
intrinsic dispersion given the large spread over two or- 
ders of magnitude (10~ 3 — 10 _1 ) in the relation between 
accret ion rate and SFR (Fig. \l3b). iRobertson et al.1 
(2006) point out that AGN samples at high redshift may 
have larger intrinsic scatter in the velocity dispersion 
of their hosts due to their immature dynamical state. 
The observational situation is s o far unclear: l i mited 
AGN samples at higher redshifts (|Woo et al.ll2006l 12008ft 
show preferentially higher black hole masses relative to 
that expected by the local M — a relation, while other 
studies of low redshift AGNs, undergoing substantial 
rates of accretion (Lb i/ 'LEdd > 0.1), have black hole 
masses lower t han that of in active galaxies of equiva- 
lent h ost mass (|Ho et al.ll2008| ) or luminosity (|Kim et al.l 
2008). Notwithst anding the large observational biases 
(|Lauer et al.ll200"7t ). it is possible that the coeval growth 
of galaxies and their SMBHs may be intermittent with ei- 
ther stars or SMBHs growing somewhat faster or slower 
than each other, as indicated by our findings exempli- 
fied in Figure [13b . while still resulting in a tight black 
hole-bulge m ass relation for in active galaxies at z = as 
suggested bv lWoo et all (|2006| ). 

6. SUMMARY AND CONCLUSIONS 

We have utilized the zCOSMOS 10k catalog of galax- 
ies with reliable spectroscopic redshifts to investigate the 
properties of those that host AGN. X-ray observations 
with XMM-Newton enable us to identify 152 AGN, from 
a parent sample of 7543 zCOSMOS galaxies, that include 
those with significant obscuration and of low optical lu- 
minosity. The derived properties of the full zCOSMOS 
sample such as stellar mass, rest-frame colors, and spec- 
tral properties (i.e., emission line strength, D„4000) en- 
able us to determine the prevalence of AGN activity as a 
function of galaxies with the aforementioned character- 
istics. 
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Specifically, we measure the SFR of galaxies using the 
[OIIJA3727 line luminosity. We account for the con- 
tribution from the underlying AGN component most 
likely arising from the narrow line region by using the 
[OIII]A5007 luminosity and the typical [OII]/[OIII] ra- 
tio found f rom previous st udies of AGN and quasars 
(|Hol [20051 : Ikim et al.1 [20061 ) . The [OIII] line luminos- 
ity is measured directly from our spectra if present. For 
the subsample with [OIII] outside our observed spectral 
bandpass, we infer the [OIII] strength from the hard (2- 
10 keV) X-ray luminosity and the well known correlation 
between these two quantities. 

Overall, we find that the two main requirements for 
a galaxy to host an actively, accreting SMBH are (1) its 
stellar mass and (2) a significant amount of gas content as 
inferred by the observed levels of star formation and the 
age of their stellar populatio ns. These findings essen tially 



igse 

extend those of the SDSS (|Kauffmann et al.ll2003bl ) out 
to z ~ 1. 



We particularly draw the following conclusions: 

• We confirm with many previous studies that the 
fraction of galaxies hosting AGN rises with increas- 
ing mass with most host galaxies having M* > 
4 x 10 10 M Q . 

• The host galaxies of AGN have ongoing star for- 
mation with a broad range of rates (» 1 — 100 
Mq yr _1 ) higher than that of the overall massive 
(log M > 10.6) galaxy population and essentially 
equivalent to those forming stars (i.e., emission- line 
galaxies). The association of AGN activity and 
young stellar populations is further substantiated 
by an observed increase in the fraction of galaxies 
harboring AGN at low values of the spectral index 
Z?„(4000) and blue rest-frame color U — V. 

• The enhancement of AGN in young, star-forming 
galaxies with 1.0 < £>„(4000) < 1.4 and 0.8 < 
U — V < 1.5 lessens the evidence for AGN prefer- 
ing to reside in galaxies of intermediate colors and 
their role in galaxy evolution. We demonstrate that 
previous studies based on luminosity-selected sam- 
ples are misleading due to the inclusion of galaxies 
with low mass-to-light ratios that are less likely to 
harbor AGN given their lower masses. 

• The SFRs of AGN hosts evolves with redshift in 
an equivalent manner to the overall star-forming 
galaxy population. This essentially brings the ev- 
idence for a co-evolution scenario between accre- 
tion onto SMBHs and the star-formation history of 
galaxies to a closer physical scale (i.e., within the 
same galaxies). 

• A direct relationship between the consumption of 
gas into stars and that accreted onto SMBHs is 
weak, suggesting that additional physical complex- 
ities or varying timescales may be inherent. On av- 
erage, a co-evolution scenario for the overall pop- 
ulation is clearly evident given the constancy of 
the ratio (~ 10~ 2 ) between mass accretion rate 
onto SMBHs and SFR with redshift possibly in- 
dicative of depleting gas reservoirs from z <~ 1 to 
the present. The order-of- magnitude increase in 
this ratio compared to the locally measured value 
of Mbh /Mbuige, is consistent with an AGN lifetime 
substantially shorter than that of star formation. 
Furthermore, the significant dispersion in this ratio 
may be indicative of larger scatter in the BH-bulge 
relations at higher redshift. 

Our results put important constraints on physical mod- 
els of AGNs and their evolution. The considerable rates 
of star formation (~ 10 M Q yr _1 ) in their hosts, coupled 
with the lack of structural signs of g alaxy interactions 
(|Grogin et al.ll2005l : lGabor et al.ll2008f) and timescales for 
accretion, indicate that a "Seyfert mode" of accretion, 
driven by secular processes, is more likely than merger- 
driven models for this class of moderate-luminosity AGN. 
In particular, there is no indication of the suppression or 
truncation of star formation at levels expected from mod- 
els implementing AGN feedback. Therefore, we find no 
conclusive evidence that AGNs, undergoing this mode of 
accretion, are a key factor in the evolution of galaxies. 
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